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Abstract 

We consider an abstract functional-differential equation derived from the pressureless Euler 
system with variable coefficients that includes several systems of partial differential equations 
arising in the fluid mechanics. Using the method of convex integration we show the existence 
of infinitely many weak solutions for prescribed initial data and kinetic energy. 
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1 Introduction 

The concept of weak solution is indispensable in the mathematical theory of inviscid fluids, where 
solutions of the underlying non-linear systems of partial differential equations are known to develop 
singularities in a finite lap of time no matter how smooth the initial data might be. The weak 
solutions are being used even in the analysis of certain viscous fluids like the standard Navier-Stokes 
system, where a rigorous theory in the classical framework represents one of the major open problems 
of modern mathematics. In the absence of a sufficiently strong dissipative mechanism, solutions of 
non-linear systems of conservation laws may develop fast oscillations and/or concentrations that 
inevitably give rise to singularities of various types. As shown in the nowadays classical work of 
Tartar pR], oscillations are involved in many problems, in particular in those arising in the context 
of inviscid fluids. 

The well know deficiency of weak solutions is that they may not be uniquely determined in terms 
of the data and suitable admissibility criteria must be imposed in order to pick up the physically 
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relevant ones, cf. Dafermos [3]. Although most of the admissibility constraints are derived from 
fundamental physical principles as the Second law of thermodynamics, their efficiency in eliminating 
the nonphysical solutions is still dubious, cf. Dafermos ms- Recently, DeLellis and Szekelyhidi HH 
developed the method previously known as convex integration in the context of fluid mechanics, in 
particular for the Euler system. Among other interesting results, they show the existence of infinitely 
many solutions to the incompressible Euler system violating many of the standard admissibility 
criteria. Later, the method was adapted to the compressible case by Chiodaroli |[7j|. 

In this note, we introduce an abstract functional-differential equation that may be viewed as the 
pressureless Euler system with variable (functionally solution dependent) coefficients. We present 
an abstract version of the so-called oscillatory lemma and use it in order to show the existence 
of infinitely many solutions adapting the method of HQ • Various specific systems arising in fluid 
dynamics will be then identified as special cases of the abstract problem. 

The paper is organized as follows. In Section [2] we introduce the abstract problem and formulate 
our main result proved in the remaining part of the paper. To this end, we adapt the apparatus of 
convex integration including the concept of subsolution in Section [3j In Section [4l we present the 
oscillatory lemma and show the existence of infinitely many solutions. Several specific examples are 
discussed in Section 0 Finally, Section [H] addresses the problem of strong continuity of the weak 
solutions at the initial time. 


2 Abstract problem, main result 


The symbol R^^ N will denote the space of N x N symmetric matrices over the Euclidean space R N , 
N = 2, 3, is its subspace of those with zero trace. For two vectors v, w e A jV , we denote 


V <g) w e R^ N , [v <g) w]jj = ViVj, and v 0 w e v©w = v<g)w- A_v • wl. 


NxN 


N 


For the sake of simplicity, we suppose the physical space to be the “flat” torus 

N 


meaning, the functions of i 6 11 are (2-)periodic in R N . 


2.1 Abstract problem 

We consider the following problem: 
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Find a vector field u G C wea k([0, T]; L 2 (f2; R N )) satisfying 
(u + h[u]) O (u + h[u]) 


chu + div 2 


r|uj 

1 |u + h[u|| 2 


+ H[u] = 0, div^u = 0 in T>'((0, T ) x 12; R N ), (2.1) 


(t, x ) = e[u](i, x) for a.a. (t, x) G (0, T) x 12, 


2 r[u] 

u(0, •) = u 0 , u(T, •) = u T , 
where h[u], r[u], H[u], and e[u] are given (nonlinear) operators. 


( 2 . 2 ) 

(2.3) 


Remark 2.1 The problem Ii2.1\ - \2.3j) is seemingly overdetermined as both the initial and the end 
state are prescribed. Moreover, the associated “kinetic energy’’ is constrained by H2.2\) . Specific 
applications will be given in Section 0 

Remark 2.2 The choice 

h = 0, r = 1, H = 0, e = e(t) 

gives rise to the pressureless (incompressible) Euler system 

d t u + div x (u ® u) = 0, div x u = 0 
with the prescribed kinetic energy 

-|u| 2 = e(t) 

2 w 

studied by Chiodaroli w and DeLellis, Szekelyhidi Ill- 
Remark 2.3 Note that a “more complex’’ problem 

(u + h[u]) (8) (u + h[u]) 


d f u + div.v 


r u 


+ H[u] + V x n[u] = 0, diva-u = 0 


(2.4) 


can be converted to \2.1\ ) . \2.2\) . with 


e[u] = Z[u](t) - yn[u], 


where Z is an arbitrary spatially homogeneous function. 
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Remark 2.4 The “pressure” II in (2-4) can be incorporated in H by solving the problem 


div^Hn = V x n in Q, H n (x) G x G CL. 

We can take, for instance, the solution of the Lame system 

2 

Hn = V X U + V X U T - — div x Ul. 

As observed by Desvillettes and Villani [1J2, Section IV.I, Proposition 11], the vector field U is 
uniquely determined up to an additive constant. Of course, in order to preserve certain continuity of 
H n , more regularity of A is needed. 

The quantities h, r, H, and e are operators depending on the solution u. In order to specify their 
properties, we introduce the following definition: 

Definition 2.1 Let Q C (0, T) x be an open set such that 

\Q\ = \(0,T) x Q\. 

An operator 

b : C' weak ([0, T\-L\n- R n )) n L°°((0, T ) x O; R N ) -G C b (Q, R M ) 
is Q— continuous if: 

• b maps bounded sets in L°°((0, T) x O; R N ) on bounded sets in C b (Q,R M ); 

• b is continuous, specifically, 

b[vn\ —> 6[v] in C b (Q ; R M ) (uniformly for (t,x) G Q ) 

whenever (2.5) 

v„ —>■ v in C wea k([0, T\; L 2 (0; R N )) and weakly-(*) in L°°((0, T) x fl; R N ); 

• b is causal (non-anticipative), meaning 

v(t, •) = w (t, •) for 0 < t < r < T implies 6[v] = b[ w] in [(0, r] x f2] D Q. (2.6) 
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In this paper, we suppose 


h = h[u] : C weak ([0, T}-L\n- R n )) -G C b (Q; R N ), 
r = r[u] : C weak ([0, T]; L 2 (L>; R N )) -G C b (Q ; R), r > 0, 
e = e[u] : C weak ([0, T]; L 2 (ft; R*)) -G <7 6 (Q; R), e > 0, 

H = H[u] : C weak ([0, T\- L\Q- R N )) -G C b (Q ; 
are given Q —continuous operators for a certain open set 0. 


2.2 Subsolutions 

Before stating our main result concerning solvability of problem (I2.lfl2.3p . it is convenient to introduce 
the set of subsolutions. Let A max [A] denote the maximal eigenvalue of a matrix A G R^ y ^ N . Similarly 
to DeLcllis and Szekelyhidi HP, we introduce the set of subsolutions: 


X n = 


£ Cwe ak ([0, T]; L 2 (h2; R N )) D L°°((0, T) x 12; R N ), v(0, •) = u 0 , v(T, •) = u T , (2.8) 


•N\ 


dtv + div x F = 0, div x v = 0 in V((0, T) x hi; iW), for some F G L°°((0, T) x fl; R^^ 0 ), 


N \ 

sup — A max 

(t,x)£Q,t>r ^ 


veC(Q-R N ), weC(Q-R^), 
(v + h[v]) <g) (v + h[v]) 


r[v] 


- F + H[v] 


— e [v] <0 for any 0 < r < T 


Remark 2.5 Note that, in contrast with m the inequality 


N 


A r 


(v + fe[v])®(v + fr[v]) 
r[v] 


- F + H[v] 


< e v 


is satisfied only on the open set Q, where all quantities are continuous. Moreover, the inequality is 
strict on any open time interval (t,T) , 0 < r < T. 


2.3 Main result 

We are ready to state our main result. 
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Theorem 2.1 Let the operators h, r, H, and e given by \2. 7| ) be Q— continuous, where Q C 
[(0, T) x fi] is an open set, 

\Q\ = |(0,T) x Q\. 

In addition, suppose that r[v] > 0 and that the mapping v i —y l/r[v] is continuous in the sense 
specified in 112.5\) . Finally, assume that the set of subsolutions X 0 is non-empty and bounded in 
L°°((0, T) x Q-R n ). 

Then problem \2.1\ - HOI) admits infinitely many solutions. 


The next two sections will be devoted to the proof of Theorem 12.11 For the set of subsolutions to 
be non-empty, the energy e must be chosen large enough. For instance, taking Uo = u-p e C(h2; R N ), 
div^Uo = 0 we check easily that X 0 is non-empty, specifically u 0 E Xq, as soon as 



max 


(u 0 + h[ u 0 ]) ® (u 0 + h[ u 0 ]) 
r[ u 0 ] 


+ H[uo] 


< e[uo]. 


Recalling the purely algebraic inequality (cf. [TT] ) 


1 h 2 N x 

-'-X- < —A 

2 f “ 2 


h ® h 

r 



(2.9) 


( 2 . 10 ) 


where the equality holds only if 


we get from (12. 9 j) that 


H 


h (g) h 1 |h| 2 
r Nr 


1 | u 0 + h[u 0 ]| 2 

2 r[u 0 ] 


< e[u 0 ], 


( 2 . 11 ) 


meaning the relation (12. 2 K is violated at the initial time. This is the undesirable initial “energy jump” 
characteristic for the weak solutions obtained by the method of convex integration. A possible remedy 
for this problem will be discussed in Section |6j 


3 Convex integration 

As the set X 0 is bounded, there exists a positive constant e such that 

e[v] < e for any v G X 0 . (3.1) 
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Under the hypotheses of Theorem 12. II we may define a topological space X 0 as the closure of the space 
of subsolutions X 0 with respect to the (metrizable) topology of C' wea k([0, T]; L 2 (12; R N )). Accordingly, 
X 0 is a (non-empty) complete metric space with the distance of two functions v, w given by 


sup d[v(t,-);w(t,-)], 
te[o,T] 


where d is the metrics induced by the weak topology on bounded sets of the Hilbert space L 2 (12; R N ). 
Note that, in view of (12.10)1 . (13.11) and boundedness of all operators involved in the definition of X 0 , 
the associated fluxes F are bounded in L°°, in particular, 


d t y + div x F = 0, div x v = 0 in V{{ 0, T ) x 12; R N ), (3.2) 

for any v G X 0 , where the flux F G L°°((0, T ) x 12; R^^, o) can be obtained as a weak limit of fluxes 
in X 0 . Moreover, by convexity of the function 



max 


v + h) ( 8 ) (v + h) 

; ------ - F + H 

r 


in v and F, we get 


N x 


(v + h[v]) <g> (v + h[v]) 
r[v] 


F + iv 


< e[v] a.a. in (0,T) x 12. 


(3.3) 


Next, we introduce a countable family of functionals 


/Jvl = 


[ T f 

1 V + h V 

/ / 

~ e v 

Jl/n Jfl 

[2 r[v] 


dx : X 0 —> (—oo, 0]. 


In accordance with the hypotheses (12.71) . each I n can be seen as a lower semi-continuous functional 
on X 0 . In particular, by means of Baire’s category argument, the set 

S = n n>0 |v G X 0 v is a point of continuity of J n | 

has infinite cardinality. 

In the next section, we show that 


if v is a point of continuity of I n in X 0 , then 7 n [v] = 0. (3-4) 

In accordance with (I2.10p . (12.lip , combined with the previous observations stated in (13.2p . (13.3p . this 
implies that S consists of weak solutions to problem (l2.1H2.3p . Consequently, the proof of Theorem 
12.II reduces to showing (13.4p . 
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4 Oscillatory lemma, infinitely many solutions 

In accordance with the previous discussion, the final step in the proof of Theorem 12.11 is to show 
(13.41) . The main tool we shall use is the following variant of the oscillatory lemma (cf. De Lellis and 
Szekelyhidi [1T[ Proposition 3], Chiodaroli (7, Section 6, formula (6.9)]) proved in [13, Lemma 3.1] : 

Lemma 4.1 Let U C Rx R N , N = 2,3 be a bounded open set. Suppose that 

h G C(U] R n ), fit G C{U ; R^ n 0 ), e, f G C(U), r > 0, e<einU 

are given such that 


N 


A r , 


h <g) h 


< e in U. 


(4.1) 


Then there exist sequences 

w„ 6 R n ), G n 6 crm R^), n = 0,1, 

such that 


N \ 
— A r 


d t w n + div x G n = 0, div 2 ,w n = 0 in R N , 
(h + w„) (g) (h + w n ) 


(fi + G n 


< e in U, 


and 


w n — > 0 weakly in L 2 (U ; R N ), lim inf J -— dx (it > A(e) J — - 
for a certain A(e) > 0 depending only on the energy upper bound e. 


dx d t 


(4.2) 


(4.3) 


Remark 4.1 Note that Lemma 4-1 applies to continuous, not necessarily bounded, functions on the 
open set U. 

With Lemma 14.11 at hand, we may show the following result that contains (13.4ft as a particular 
case. 


Lemma 4.2 Let 


In = 


ID 


1 |v + h[v] 

2 r[v] 


e v 


dx dt : Xq —> (—oo, 0] 


be a functional defined on an open set D C [(r, T) x fl] n Q, 0 < r < T. 
Then Id vanishes at any of its points of continuity. 






















Proof: 

Arguing by contradiction we assume that v e X 0 is a point of continuity of Id such that 

J D [v] < 0. 

Since Id is continuous at v, there is a sequence {v m }“ =1 C X 0 (with the associated fluxes F m ) such 
that 


v m ->• v in (7 weak ([0,7^;L 2 (n;i2 JV )), I D [v m ] Id[v] as m oo. 
As [v m ,F m ] are subsolutions and r > 0, we get, thanks to 


-A 

2 ^max 


< elv r 


h[v m ]) ® (v m + h[v m ]) 
r[vj 


— F m + M['' 


■ 8 m in [[r, T) x 0)] D Q for some 8 m \ 0. 

Now, fixing m for a while, we apply Lemma [4.II with 

N = 2,3, U = D, r = r[v m \, h = v m + h[v m ], H = F m - H[v m ], e = e[v m ] - S m . 

For {[w min , G min ]}^ =1 the quantities resulting from the conclusion of Lemma 14711 we set 


Vm,n "L W mjn , F }7ljn F m Gtyn^n- 


Obviously, 

‘dtym,.n A div x F mjn 0, di VxV m ,n 0 in ID ((0, F) x O), v mri (0, ■) Vo, v miri (7~', •) vy. 
Moreover, in accordance with (14.2(1 and the fact that w n , G m , n vanish outside D, 


N 


Xr 


f m,n 


h[v m ]) (8) (v m>n + h[v m ]) 


r W 


- F m>n + H[v m ] 


< e[v m ] - 5 m in [r, T) x O 0 Q, 


and, by virtue of the causality property (12.611 . 


N \ 

sup — A max 

(t,x)(E:Q,S<t<T ^ 


h[v m?n ]) (v m,n h["V m?n ]) 


F m ,n + 




for any 0 < s < r. Consequently, in view of continuity of the operators v i— y h[v], r[v], e[v], H[v] 
specihed in (12.5(1 . we may infer that for each m there exists n = n(m) such that 


^m,n(m) ^ Xq, 171 1, 2, 
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Moreover, by virtue of (14.3p . we may suppose that 


v m ,„( m ) -t v in C' weak ([0,T];L 2 (O; R 2 )) 


in particular, 

lD[v m ,n(m)] ^ Id[v] (4.4) 

as m —> oo. 

Finally, using again the conclusion of Lemma [4.11 combined with Jensen’s inequality, we observe 
that the sequence can be taken in such a way that 


lim inf / 

m—> oo 


D m,n(rn 


= lim inf 

' J m— 


D 


( 1 | 4“ W m,n(m) 4" h[v m 4“ W m,n(m)\ \ 

\2 r[v m 4 - W mitl(m) ] 


e[v m 4- w TO) „ 



dx dt 


/1 |v m 4- h[v m 4- w m;n(r] 

™ 4 °° Jd V 2 r [V m + W m ,n{m) ] 


= lim 


e[v m + w m>n ( m )] dx dt 


lim inf 


1 |w 


- T dx dt 


Jd 2 r[v m 4- W min(m) ] 


> I D [w} + lim inf f ( e([v m ]) — S m — 

2 m—Hx J D y 

- /d|v) + 1 iSS f ( J D - S ’ 1 |Vm + hk 


1 |v m + hjv^ 

2 r[v m ] 

112 \ \2 


12 \ 2 


dx dt 


2 r[v m ] 


which is compatible with (14. 4 ft only if Id [v] = 0. 


dx dtj = I D [v] + (Jd[v]) 2 , 


Q.E.D. 


We have shown (13.4j) : whence Theorem 12.11 


5 Examples 

There are many systems arising in mathematical fluid dynamics that can be written in the abstract 
form (12. II [273]) . We review some of them already studied in the available literature. 
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5.1 Euler-Fourier system 


The Euler-Fourier system describes the time evolutions of the mass density g = g(t,x), the velocity 
u = u (t,x), and the (absolute) temperature $ = $(£, x): 

d t g + div x (£>u) = 0, (5.1) 

d t (gu) + div x (^u <g> u) + V x (gtf) = 0, (5.2) 

^(d t (gi9) + div x (^u)) - Ai9 = -gd div x ii. (5.3) 

Following [ 8 ] we first write the momentum £>u as its Helmholtz decomposition 

0 u = v + V^, div x v = 0 . 

Accordingly, we may fix the density g and the acoustic potential $ so that 

d t g + A<f> = 0 holds, 

meaning equation (15.11) is satisfies as div x v = 0. 

With g, <3> given we may determine the temperature field $ = $[v] as the (unique solution) of 
(15.31) , specifically 


{gd t d + (v + V X T) • V x $) — A$ = — 0 iMiv x (v + 


V x 4>) 


endowed with appropriate initial data. 
Finally, we rewrite (I5.2j) in the form 


div + div x | ( v + V ^)^( v + V ^) ^\ + Vx («%V X $ + g$[v}) = 0. 


Fixing the “energy” so that 

1 |v + V x $| 2 


Q 


N 

= e[v] = Z - — (dtV x $ - [v]), 


2 r 

where Z — Z(t) is a suitable spatially homogeneous function, we reduce (I5.4[) to 

d t v + div, ( (v + V.«)8(v + V,«) \ = 0 _ d . viv = 


(5.4) 


(5.5) 


(5.6) 


which is an equation in the form (12.ip . 

With certain effort, it is possible to show that the hypotheses of Theorem 12.11 are satisfied for 
Q = (0, T) x O, and we obtain the following result, see [81 Theorem 3.1]: 


11 












Theorem 5.1 Let T > 0 be given, along with the initial data 


0 ( 0 , ■) = 00 e ft, > 0, 0(0, •) = 00 e C 2 (fi), 00 > 0, u(0, •) = u 0 e C 3 (fi; R N ), (5.7) 

H = ([-1;1] { _ 1;1} )", N = 2,3. 

Then the Euler-Fourier system < 15.11 - 15.51) admits infinitely many weak solutions in (0, T) x 
emanating from the same initial data d 5. 7|). 


As already pointed out, the solutions obtained in Theorem 15.11 mav be non-physical in the sense 
they violate the principle of energy conservation. However, this drawback can be removed at least 
for certain initial data. We will discuss this issue in Section 16.11 


5.2 Quantum fluids 

The Euler-Korteweg-Poisson system describes the time evolution of the density g = g{t, x) and the 
momentum J = J (t,x) of an inviscid fluid: 


d t J + dh^ ^ J * J j 


+ V x p(ft) 


d t g + div x J = 0, 

-aJ + gS7 x (jC(g)A x g + ^K'(g)\W x g\ 2 ^j + qSJ x V, 

= 0 - ft 


(5.8) 

(5.9) 
(5.10) 


where K : (0, oo) —> (0, oo) is a given function, see Audiard |3] , Benzoni-Gavage et al. [4] , [5j. The 
choice K = K > 0 yields the standard equations of an inviscid capillary fluid (see Bresch et al. 
[6], Kotchote [16], [T7j), while K(g) = 4 gives rise to the quantum fluid system (see for instance 
Antonelli and Marcati m, ra. Jiingel [151 Chapter 14] and the references therein). 

For 

X(g) = gK(g ), 


it can be shown that system (15.81 [5TT0h can be recast in the form 


d t v + div 2 


(v + h) 0 (v + h) 


+ h + v x n = o, 


(5.11) 


with 


r = e f g, h = e f V x M, 

1 , 1. 


H(M) = 4e* (x(q)^xVq® gX(0)|V a .v / 0| 2 I - ® V X V + —\\7 x V\ 2 Ij , 


12 















and 


n (t,x) = e t (p(g) +d t M + M - x(q)A x Q - ^x'(q) I V x p| 2 + t^x(q) |V xx /^| 2 - gV + ^|V x y| 2 ^ , 

where g and M are suitably chosen functions, see [IB] . 

Now, Theorem 12.11 can be applied to obtain the following result, see [13, Theorem 2.1] and the 
proof therein. 


Theorem 5.2 Let T > 0 be given. Suppose that p and x satisfy 

p G C 1 ^, oo) D C 2 (0, oo), p(0) = 0, x G C 2 [0, oo), x > 0 in (0, oo). 

Let the initial data be given such that 

p(0, •) = Qo = r'o, r 0 G C 2 (h2), rneas jx G r 0 (x) = o| = 0, (5-12) 

J( 0 , •) = J 0 = p 0 Uo, U 0 G R 3 ). (5.13) 

Then the initial value problem \5.8W5.1(K) . L5.12\) . H5.13 1) admits infinitely many weak solutions 
in (0, T ) x Q. 

In the situation described in Theorem 15.21 the set Q must be taken 

Q = (0,T) x Q\{(t,x) £>(t,x) = oj. 


5.3 Binary mixtures of compressible fluids 

We consider a physically motivated regularization of the Euler equations proposed in the seminal 
paper by Lowengrub and Truskinovsky |18j . The model describes the motion of a mixture of two 
immiscible compressible fluids in terms of the density g = g(t,x), the velocity u = u (t,x), and the 
concentration difference c = c(t,x). The fluid is described by means of the standard Euler system 
coupled with the Cahn-Hilliard equation describing the evolution of c: 


d t g + div x (gu) = 0 , 

d t (gu) + div x (pu <8) u) + V x p 0 (p, c) = div x ^ gS/ x c ® V x c - || V x c| 2 I^ , 
d t (gc) + div x (pcu) = A (p 0 (Q, c) - ^div x (qX7 x c) j , 


(5.14) 

(5.15) 

(5.16) 
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where 


(5.17) 


Po(q,c) = g 


2 d fo{Q,c) df 0 (g,c) 


" i ho (Qi C) 


dg dc 

for a given free energy function /q. The system is neither purely hyperbolic nor parabolic as the 
dissipation mechanism acts in a very subtle way through the coupling of the Euler and the Cahn- 
Hilliard systems. 

The machinery of convex integration can be applied, first fixing g and <E>, similarly to Section [5711 
to solve 

d t o + A<f> = 0, 

then taking c = c[v], div x v = 0 to be the unique solution of the equation 


d t (gc) + div a; (^c(v +V x <h)) = A f/r 0 (h,c) - -div x (gV x c) J . 


Accordingly, we obtain 

d t v + div^, 
1 |v +V x $| 


v(0, •) = v 0 , div x v = 0, 
(v + V x $) O (v + V x $) 


6 


= e[v] = Z(t) 


N (l 


2 \6 


- Q (Va,.c[v] © V x c[v]) = 0, 


I V ; cc[v]| 2 + po(g, c[v]) + <%v x <h 


where Z is a spatially homogeneous function. 

Theorem 12 .II yields the following result, see jT4j for details: 

Theorem 5.3 Let the potential /o = fo(g, c) satisfy 


1 — C 1 -J- c 

fo(g, c) = H{c) + log(^) ( ai—-h «2 


(5.18) 

(5.19) 

(5.20) 


. _ , H e C 2 (R), \H"(c)\< H for allce R 1 . 

Then for any choice of initial conditions 

0(0, •) = ho e C 3 (fi), inf p 0 > 0, u(0, •) = u 0 G C 3 (fi; A 3 ), c(0, ■) = c 0 <E C 2 (fl), 


the problem d 5.1f\ - 15 . 1 61) admits infinitely many weak solutions in (0 ,T) x 
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6 Continuity at the initial time, admissible solutions 

The major drawback of the construction delineated in the previous part of the paper and the main 
reason why the weak solution obtained via convex integration can be eliminated as physically unac¬ 
ceptable is the energy jump at the initial time discussed in Section I2T31 On the other hand, however, 
once a subsolution v along with the associated energy e[v] are obtained, it is possible to show the 
existence of another subsolution defined on a possibly shorter time interval for which the initial en¬ 
ergy is attained. Such a subsolution can be then used in the process of convex integration to produce 
weak solutions that are strongly continuous at the initial time and dissipate energy. 

We first state the result for the abstract system and then shortly comment on possible applications. 
Modifying slightly the procedure used in the proof of Theorem 12.II we can show the following assertion: 


Theorem 6.1 In addition to the hypotheses of Theorem \2.1l suppose that 


|a: G Q (t,x) G <5}| = |f2| for any 0 < t < T. 


( 6 . 1 ) 


Then there exists a set of times 7Z C (0,T) dense in (0,T) such that for any t E IZ there is 
v G X 0 with the following properties: 


G a([(0, r) U (r, T) x Q] n Q; R N ) n C' weak ([0, T]-,L 2 (Tl; R N )), v(0, •) = u 0 , v(T, •) = u T ; 


2 /n. dAG 


( 6 . 2 ) 


dt'v + div x F = 0, div x v = 0 in ^((O, T ) x fl; R N ) 


(6.3) 


for some F G C b ([(0, r) U (r, T) x fl] fl Q; i^ym, 0 ); 


N 


-A, 


(v + fe[v]) ( 8 ) (v + fe[v]) 
r[v] 


F + H v 


< e[v] in [(0, r) x 12] fl Q, (6.4) 


sup — A max 

(t,i)eQ,i>r+s 2 


(v + fe[v])<8)(v + fe[v]) 
r[v] 


F + H[v] 


e[v] <0 for any 0 < s < T — r, 


(6.5) 


7 J V+ Hvl V " G,-)ds = / n e[v](T,.)dr. 


2 Jn r[v] 


( 6 . 6 ) 
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Remark 6.1 Unlike the subsolutions considered in the proof of Theorem \2.1[ the function v satisfies 
46'. 61) and is therefore strongly continuous at the point r attaining the desired energy e[v](r, •) in the 
integral sense. There is no energy jump at the time t — t! Moreover, the set of such times is dense 
in (0, T). 

Remark 6.2 In view of hypothesis m we have 

C b (Q ; R m ) C Cioc(0, T; L q (tt: R M )) for any 1 < q < oo. (6.7) 

This observation justifies 46.61) and will be frequently used in the proof below. 

Proof: 

The function v will be constructed recursively as a limit 

~v k —y v in C wea k([0, T\] L 2 (f2; R N )) for suitable v fc G X 0 . (6.8) 

We start by fixing the open interval (ao, 6o) C (0, T) in which the time r is to be localized. As 
the space X 0 of subsolutions is non-empty, we take 

v 0 € Ao, along with the associated flux F 0 . 

Next, we construct a sequence of functions v*,, open intervals (a*,, bk) C (0, T), times Tk G (a*,, bk), 
and a decreasing sequence of positive numbers Sk \ 0 such that: 


d t v k + div x F fc = 0, div x v fc = 0 in V'((0,T) x f2), v(0) = u 0 , v(T) = u T , 
for a certain field ¥ k G C(Q\ Rf y 


v fc - v fc _i G Cf°(Q] R N ), supp[v k - v fc _i] c [(a k , b k ) xfiJflQ, 
F fc - F fc „i G C£°(Q; R^ N 0 ), supp[F fc - F^] C [(a k , b k ) X fi] n Q, 

where 

0 < a k - 1 < a k < b k < b k - 1, e k = b k - a fc ->■ 0 for k ->■ oo; 


(6.9) 


( 6 . 10 ) 


sup d(v*(t), v fc _!(t)) < (6.11) 

te[o,T] 2'- 
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sup 

te(o,T) 


< 


2 k 


[ -r—r( v fc - v fc _i) • v, dx 

r[ V j] 

for all j — 0 ,..., k — 1 ; 

there exists G (a*, 6 *,) and a positive constant A independent of k such that 


1 r jvfc + h[v fc ]|- 

2 Jn r[v fc ] 

1 f |v fc _i + h[v fc _i]| 2 


(r fc , ■) dx>- [ 

2 Jf! 


1 f |v fc _i + h[v fc _ 


o r[v fc _i] 

1, 1 




> x / ' “ " r — L r — (r k - i, •) dx + ~X—a k for all t G (a fc , 6 fe ), 
^ jo r’|v/ c _;ij 2 £ k 


where we have set 


Oik — 


I L f e t vfc - 

J CL fc «/ Q V 


1 v fc _i + h v fe _i 2 A 

11 - o-r-i- dxd t > 0; 

2 r Vfe-i / 


( 6 . 12 ) 


-— (t, •) dx + A-j(6.13) 


iV, 


TV 


An 


(Vfc + h[y fc ]) ( 8 ) (Vfc + h[v fc ]) _ 

r[v k ] 

(Vfc + h[vfc]) (8) (Vfc + h[vfc]) 


Ffc + H[vfc] 


’[Vfc] 


- Ffc + mjvfc] 


< e[vfc] in [(0, a fc ] xfiJnQ, (6.14) 

1 


< e[vfc] — <5fc ^1 + —^ (6.15) 


in [(a fc , fefc) x 0] fl Q, 


N 


A, 


(Vfc + h[vfc]) (8) (vfc + h[vfc]) 
r[vfc] 


- Ffc + H[vfc] 


< e[vfc] - Sj_ i (l + 


m 


[[bj, bj_ i) x fl] fl <5, j = 0,1,..., k, 6_1 = T. 


(6.16) 


Step 1 

It follows from the properties of X 0 that v 0 satisfies (16.91) . along with the bounds (16.141 16.161) 
for a certain 


<5o = 5-i > 0. 


Step 2 

Suppose we have already constructed the functions Vj, along with intervals ( dj,bj ), the times Tj, 
and the constants Sj, for j — 0 , 1 ,... k — 1 enjoying the properties (16.91 - 16961) . Our goal is to find 
Vfc, (o/c,5fc), Tfc, and 5fc. 
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First, we fix the interval ( a k ,b k ). To this end, compute 

= C l ( e|Vt - l] “ \ 

As a consequence of (16.71) . the integrand is a continuous function of time continuous function of time; 
whence 

- = - f k / [eW-A ~ \ 

£k ^ a k ^Q \ 2 

f ( r i 1 |v fc _i + h[v fc _!]| 2 ^ 

J ^ I e[v fc _ij - -- ^ -I dx as e k = b k - a k ->• 0. 

Consequently, keeping in mind that a*, > 0 and repeating the same continuity argument, we may 
choose a k ~i < a k < b k < b k ~ i and e k so small that 


jvfc-i + h[vfe-i 

r[v k -i] 


dxdt 


I Vfc— i + hjyfc-i 

r[v k -i] 


dxdf. 


J_ [ bk f + h[v fe _r] 

fa J CLfc J 2 


> 


r[v k - J 
1 |v fc _i + h[v fc _i]| 2 


dxdt + A(e) 




In 2 


r[v fc _! 


> 


1 |vfe_i + h[vfc_i 


(V) dx +A(e)^ 


(r fc _i, •) dx + A(e)—j for all t e (a k , b k ), 


K 


Jn 2 rfvfc-i] 

where A(e) is the universal constant introduced in Lemma [4.11 
At this stage, we apply Lemma [4. II for 


(6.17) 


U = [(a k , b k ) x O] fl Q, h = v fc _i + h[v fc _i], r = r[v fe _i], H = F fc _i - H[v fc _i], 


and 

e = e[v fc _!] -8 k (l + ^zrj , 

where S k > 0 is chosen small enough so that (14.ip may hold. 
Now we claim that it is possible to take 


v fc = v fc _i + w„, F fc = F fc _i + G n , n large enough, 


(6.18) 


where w n , are the quantities constructed in Lemma [4.11 Obviously, the functions v k satisfy (16.91 
- 16.121) provided n is large enough. Indeed we observe that n can be chosen so large for (16.121) to be 
satisfied. To see this we realize that, by virtue of (16.7|) . the image 


U te[a k ,b k ] 


■(t, •) is compact in L 2 (0; R N ), j — 0,..., k — 1. 


r Vo- 
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Next, we use continuity of the operators h, r specified in (16.71) to compute 


r>, f 1 K + h[v,r dx d( = [<■* 


+ 


' a k 


rbk 


In 2 


r[v k ] 
2 


1 Vfc—! + h[v*_, + wjp dx dt 


/ a/g 


w, 


dx dt + 2 


rH 


In 2 r[v fc _ 1 + w n ] 

/' W n • (v fe _i + h[v fc _i + w n ]) 


a k Jn 2 r[v fc _! + w„] 

r bk r 1 jy fc _i + h[v fc _ 112 

Jfi 2 




r[v fc _! + w n ] 


dx dt 


> 


lj 1 dx dt + 7 - 


rbk 


1 |w, n 


dx dt + 


r[v fc _i] 2Ja k Jn2r[v k _ 1 \ 

provided n is large enough, where e n —> 0 as n —> oo for any hxed k. Consequently, (14.3p gives rise 


to 


/ <2fc 


“ f ilniTM! ix di + e „ > O, _ 

f 2 2 r[vfc] Ja k Jn 2 


1 |v fc _i + h[v fc _i]| 2 


r[vk-1 


dx dt 


(6.19) 


A(e) f b k 


r l ( e|vt - i] - 2 


1 |v fc _i + h[v fc _i 


r[v k - 1 ] 


dx dt 


> 


rb k 


/ <2fc 


1 |v fc _i + h[v fc _! 


In 2 




dx dt 


A(e) 1 a\ 
4 |0| e-fc ’ 


where the last line follows from Jensen’s inequality. Thus, using (I6.17I) , (I6.19p . we may find n large 
enough and r k E (a k , b k ) such that (16.131) holds with some A that can be determined in terms of A(e) 
and |fi|. 

Finally, our goal is to check that v*., satisfy (I6.14l - l6.16p . First we claim that (16.14p is a direct 
consequence of the causality property (12.61) . Next, Lemma (14.ip . specifically (14.2)) . yields 



max 


(Vfc + h[v fc _i]) ® (Vfc + hjvfc-i]) 
r[v k - 1 ] 


Ffc + H[v fc _i] 


< e[v fc _i] - S k 



in [(a k ,b k ) x hi] D Q] whence (16.151) follows from uniform continuity of h, r, H and e provided n is 
large enough. To see (I6.16p . we have to realize that v k = v k _i and ¥ k = F*,_i in [(b k ,T) x fi] D Q, 
and, similarly to the above, relation (16.16p follows from continuity of h, r, H and e as son as n is 
chosen large enough. 

Step 3 

Our ultimate goal is to observe that v, determined by the limit (16.8p . enjoys the desired properties 
claimed in Theorem 16.11 We set r = lim*.^,*, Since the functions v k , ¥ k coincide with Vfc_ 1; Ffc_! 
on the time intervals (0,a.*), ( b k ,T ), the properties (16.21 - 16.5p follow by taking the limit in (16.91) . 
(16.141 I6.16P for k oo. 
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To see (16.61) . we first observe that, by virtue of (16.131) . 


I Vfc + h[v fc ]p 


{jk) dx /* Y as k —> oo, 


2 Jn r[v fc ] 

and the convergence is uniform on the time intervals (a*,,^); whence 


_ J_ f bk 


I L ~ 


1 jVfc_L + h[v fc _! 

2 r[v fc _i] 


dxdt —> 0, 


(t) dx —> / e[v](r) dx as k —» cx> uniformly for t G (a*,, &*,). 
Jo. 


which in turn implies 

1 r |v fc + h[v fc ]| 2 

2 r[v fc ] 

We show that (j6.2Q[) yields 

w fc (r, •) ->■ w(r, •) 

which completes the proof of Theorem 16.11 Indeed we may write 


( 6 . 20 ) 


f |Vm ~ Vn| 

In r[v n ] 


dx = 


dx 


iv r ,r . _ r (v m - v„) • v n 


o r[v n ] Jn r[v n ] 


dx — 2 


'a r[v n ] 


dx, m > n, 


where the difference of the first two integrals vanishes as n —» oo uniformly for t G (0,T); whereas 


V n ) * V n 


r[v n ] 


dx 


y [ hui: dI 




r[v n ] 


< )) _ 1 uniformly in (0, T) 


in view of ( 16T2D . 


Q.E.D. 


Now, we can define a set of subsolutions on the time interval (r, T), with 

u 0 = v(t), Q = Qt = Q n [(r, T) x Q ], 
and the operators h r , r T , e T , H r defined as 

v in [0, r] 


h T [w] = h[w]|( T) T), where w = 


w G [r, T} 
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where v is the function constructed in Theorem 16.11 In accordance with (16. hi) , we have vL^] is a 
subsolution, and 


If 

2 Jn 


I u 0 + h r [u 0 ] p 
r T [ u 0 ] 


dx = / e T [u 0 ] dx. 
Jn 


Finally we note that in this case the weak solutions u constructed via Theorem 12.11 will satisfy 


1 |u + h[u]| 2 

2 r[u] 


(t,x) 


e[u](t,x) for a.a. ( t,x ) G (0,T) x 0 and including the initial time t = 0. 


6.1 Example, dissipative solutions to the Euler-Fourier system 

Revisiting the Euler-Fourier system introduced in Section 15.11 we say that g, d, u is a dissipative 
solution of (15.11 - I5T>|) . if, in addition, the energy balance 

E(t) = Qp|u | 2 + ^gd^j (t, •) dx = (^o|u 0 | 2 + l g 0 A) dx (6.21) 

holds for a.a. t G (0,T). 

As a possible application of Theorem 16.11 one can show the following result, see [8] Theorem 4.2]: 
Theorem 6.2 Under the hypotheses of Theorem \5.1l let T > 0 and the data 

go, $o £ go,‘do > 0 


be given. 

Then there exists u 0 G R N ) such that the Euler-Fourier system \5.1\ - 15.31) . with the initial 

conditions \5. 7|) ; admits infinitely many dissipative solutions in (0, T) x fh 
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